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O ■ Abstract 

"*^h Hanle scattering is an important diagnostic tool to study weak solar magnetic fields. Partial ffe- 

Q\ . quency redistribution (PRD) is necessary to interpret the linear polarization observed in strong reso- 

nance lines. Usually angle-averaged PRD functions are used to analyze linear polarization. However 
it is established that angle- dependent PRD functions are often necessary to interpret polarization 
f^> ' profiles formed in the presence of weak magnetic fields. Our aim is to present an efficient decom- 






position technique, and the numerical method to solve the concerned angle-dependent line transfer 
problem. Together with the standard Stokes decomposition technique we employ Fourier expansion 
over the outgoing azimuth angle to express in a more convenient form, the angle-dependent PRD 
function for the Hanle effect. It allows the use of angle-dependent frequency domains of Bommier 
to solve the Hanle transfer problem. Such an approach is self-consistent and accurate compared 
to a recent approach where angle-averaged frequency domains were used to solve the same prob- 
lem. We show that it is necessary to incorporate angle-dependent frequency domains instead of 
angle-averaged frequency domains to solve the Hanle transfer problem accurately, especially for the 
Stokes U parameter. The importance of using angle-dependent domains has been highlighted by 
taking the example of Hanle effect in the case of line transfer with vertical magnetic fields in a 
slab atmosphere. We have also studied the case of polarized line formation when micro-turbulent 

£N) ' magnetic fields are present. The difference between angle-averaged and angle-dependent solutions 

C*") , is enhanced by the presence of micro-turbulent fields. 

in 



o 



a 



Keywords: Line: formation - polarization - scattering - magnetic fields - methods: numerical 



1. Introduction 



The polarization of line radiation is caused by resonance scattering on bound atomic levels. 
A modification of this process by external magnetic fields is called the Hanle effect. The linear 
polarization of strong resonance lines, are particularly sensitive to the type of the frequency re- 
distribution mechanism used in their evaluation especially in the presence of magnetic fields. The 
differences of the diffuse radiation field between the linear polarization (Q/I) profiles computed 
using angle-averaged and angle-dependent partial frequency redistribution (PRD) functions are il- 
lustrated in Faurobert [1| and Sampoorna et al. [2[ in the non-magnetic (Rayleigh) case. Nagendra 
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et al. [3] showed that Stokes U profiles computed in planar slabs, for the case of Hanle effect, using 
the angle-averaged PRD functions differ significantly from those computed using angle-dependent 
PRD functions. 

In the case of angle-dependent PRD functions, the strong coupling that exist between the in- 
coming and scattered radiation makes their evaluation and subsequent use in transfer equation 
numerically expensive. The use of decomposition technique developed by Frisch y| for the Hanle 
effect and Frisch [5| for the Rayleigh case, simplifies this numerically expensive problem. In the 
non-magnetic case Sampoorna et al. |2J used this decomposition technique and developed numeri- 
cal methods to solve the polarized transfer problem with angle-dependent PRD functions . They 
also present a detailed historical account of the works on angle- dependent PRD in spectral line 
polarization. In Sampoorna [6[ Hanle transfer problem with angle- dependent PRD was solved us- 
ing single scattering approximation. Further in Nagendra & Sampoorna [7| (hereafter, NS11) the 
full Hanle transfer problem with angle-dependent PRD functions was solved by including multiple 
scattering terms and using scattering expansion method (SEM). It may be noted that SEM was 
first formulated by Frisch et al. [8| for solving the polarized line transfer equation with complete 
frequency redistribution (CRD). 

In all the above mentioned papers, a Fourier-expansion of the angle-dependent PRD function 
over azimuth angle difference (x — x') is employed, where x an d x' are the azimuth angles of the 
outgoing and incoming rays. This technique was first introduced by Domke & Hubeny [9| and was 
further developed by Frisch [J] for the Hanle transfer problem. The decomposition technique of 
Frisch 4] allowed the Hanle transfer problem to be solved in an azimuth independent Fourier basis. 
In NS11 this decomposition technique was used together with the angle- averaged frequency domains 
(approximation III of Bommier |lOj) to solve the Hanle transfer problem with angle-dependent PRD 
functions. As one has to use in principle, angle-dependent domains themselves for angle-dependent 
PRD transfer problems, the approach taken in NS11 is inconsistent. Such an approximate ap- 
proach was taken, because, it allowed to work in an azimuth independent Fourier basis. Clearly 
in an axisymmetric Fourier basis one cannot apply angle-dependent frequency domains as they 
explicitly depend on azimuth difference (x — x')- This inconsistency is at the base of the slight 
differences in the angle-dependent Hanle solutions presented in NS11 (see Fig. 4 in that paper), and 



those presented in Nagendra et al. [3J. Indeed Anusha & Nagendra [Tl| pointed out that the use 
of angle- averaged frequency domains for the angle- dependent Hanle transfer problem (as done in 
NS11) results in a loss of information. 

To overcome this inconsistency we adopt Fourier expansion of angle-dependent Hanle PRD ma- 



trix over only the outgoing azimuth angle x as suggested in Anusha & Nagendra [ll| (see also 
Anusha & Nagendra [12J). Such an expansion was proposed to solve polarized transfer problems 
in multi-dimensional media, where the radiation field is non-axisymmetric even in the absence of a 
magnetic field. In this paper we apply the decomposition proposed by them to the simpler case of 
polarized transfer in ID media, and in the presence of a magnetic field. We show that an expansion 
only over \ allows to incorporate angle-dependent frequency domains for angle-dependent PRD 
functions 'self consistently' to solve the transfer problem in the magnetic case. 

In Sect. 2 of the paper we describe the decomposition technique employed for the Hanle effect 
with angle-dependent PRD. In Sect. 3, we discuss the behavior of azimuthal Fourier components 
of redistribution matrix elements. In Sect. 4, we give the equations of the SEM to solve the Hanle 
transfer problem. In Sect. 5, we discuss the results obtained by considering our new method of az- 
imuth expansion. A comparison with the results obtained from the perturbation method described 
in Nagendra et al. [3j] and those obtained by using angle- averaged domains (in NS11) is done. In 



the same section we also revisit the well known problem of vertical field Hanle effect which arises 
only due to the angle-dependent PRD in line scattering. Further, we discuss in detail the role of 
micro-turbulent magnetic fields on line transfer using angle-averaged (AA) and angle-dependent 
(AD) versions of the redistribution matrix. Conclusions are presented in Sect. 6. 

2. The decomposition technique 

The polarized transfer equation for the Stokes vector can be written in the component form as 

31 
»-£ = \p(x) + r] [Ii( T , x, O) - Si(r, x, O)] , (1) 

where i = 0, 1, 2 refer to the Stokes parameters (/, Q, U) respectively. The ray direction is given by 
n = (0, x), with 6 — cos _1 (/j,) and \ being the polar angles, x is the frequency in non-dimensional 
units. The line optical depth is denoted by r and <p(x) is the normalized Voigt function H(a, x), 
where a represents a constant damping parameter. The ratio of continuum to the line absorption 
coefficient is denoted by r. The total source vector is given by 

Si(TiI) n ) = ^M^Mi ) (2) 

(f{x) + r 

where S c ,i are the components of the unpolarized continuum source vector. We assume that S c ,o — 
B VQ , where B Vo is the Planck function at the line center, and S c ,i = S c> 2 = 0. The line source 
vector can be written as 



5t,.(r,x,n) - C ,(r) +rt± tei^S J,(r.xMr)g^, 



(3) 



where il' (9 1 , \') is the direction of the incoming ray defined with respect to the atmospheric normal. 
The solid angle element dfi' = sin 6' A9' d%' where 8' G [0, 7r] and %' G [0, 2tt]. The primary source is 
assumed to be unpolarized, so that Go(r) = eB Vo and G\{t) = C?2(t) = 0. Here, Rij(x, tt, x',il', B) 
is the Hanle redistribution matrix with angle-dependent PRD, and B represents an oriented vector 
magnetic field. The thermalization parameter e = r//(r^ + r/), with Tj and Tr being the inelastic 
collisional de-excitation rate and the radiative de-excitation rate respectively. 

In the decomposition method used in this paper, the Stokes vector (/, Q, U) is first decomposed 
into a set of six irreducible components Iq , using which we can construct an infinite set of integral 
equations for their Fourier coefficients. Following Frisch [4| we can decompose the Stokes source 
vector into six irreducible components Sq as 

__ Q=+i< 

s i (r,x,n)= Y, E T£(i t n)sg(r t x t n), (4) 

#-=0,2 Q=-K 

with a similar decomposition for the Stokes vector Ii in terms of Iq. The Tjf(i, fi) are irreducible 



spherical tensors for polarimetry introduced by Landi Degl'Innocenti 13]. The irreducible line 
source vector components are then given by 



5f Q (r, x, fi) = G«(t) + J + " dx' j *£L 



Q"=+K' 

x^E E TlQ$(m,x,x',Q,B)I^(r,x',n'), (5) 

^ ' K'=0,2 Q" = ~K' 

where 6 [= cos _1 (n.n')] is the scattering angle, G^(t) = (eB„ , 0,0, 0,0, 0) T , and 7^|",' are the 
elements of angle-dependent Hanle redistribution matrix given by 

Q'=+K 

H%§,{m,x,x',Q,B)= Y, W$ Q ,{m,B)R n {x,x',Q) 

Q'=-K 

+ Af% Q ,(m,B)R m (x,x',Q)}T KQ , !K , Q „(n'). (6) 

In this paper we use approximation II of Bommier [l0( according to which the redistribution matrix 
is written as a product of magnetic kernel Mqq, (m, B) and the angle-dependent redistribution 
functions Rii t m{x, x' , 0) of Hummer [14J. Here the index m (= 1, 2, 3, 4, 5) stands for different 
frequency domains which depend on (x, x' ', 0). The coefficients Txq',k'Q"(^') arc defined by 

The irreducible components Iq and S'q and the coefficients Tkq',K'Q" are complex quantities. For 
practical computations, we prefer working with the real quantities. In order to transfer complex 
quantities into the real space, we follow the procedure given in Frisch |l5| . First we define 

I«> K ( T ,x,n) = Re{I%(T,x,n)}, 

I$> y (T,x,n) = Im{I$(T,x,n)}. (8) 

Using these real components it can be shown that <S r = (S®, Sq, S' 1 ' x , <Sy y , S 2 ' x , S 2 ' y ) T and the 
corresponding intensity vector I r satisfy the transfer equation given in Eq. (flj with Si and 1% 
replaced by <S' and X r respectively. Now using the T matrix given in Section 5.3 of Frisch [15j the 
irreducible line source vector in terms of the real quantities can be written as 



/+oo P 

dx'f 



47T 

x "iy E E Tl^'(m,x,x',e,B)I^'(r,x',il'). (9) 

Here 7Z-qq„ has the same form as Eq. ([6]) with Nqq, and Tkq',k'Q" replaced by Nqqi and 
T^-g, k iqii. The Q indices take values [0, +K]. The elements of matrix T y K q, K >Qt>(fi') are listed 

in Appendix D of Anusha & Nagendra [16J. The explicit form of -Mqqi (m, B) can be found in 
Appendix A of Anusha et al. [17] where they are denoted by M^ l \B) with % playing the role of m 
in our notation. The formal solution of the transfer equation can now be written as 

>0, 



% K (j,x,n)= f DC e - (r '- r) ^ )/ ^ Q ^(T,x,n)^dr', for/i 
Jt ' A* 

I 1 A K (t, x,tt) = - I' e-V-T)<P<.*)li* S"a K {t, x, «)^ dr', for ,i < 
w Jo v A* 



(10) 



The irreducible components of the line source vector in Eq. ([§]) continue to be non-axisymmetric, 
because of the presence of angle-dependent redistribution function. It is computationally advan- 
tageous to express S^'q in terms of axisymmetric irreducible components. This can be achieved 
through the introduction of Fourier azimuthal expansion of the angle- dependent PRD functions. 
In this paper we use approximation II of Bommier [lOj , the expressions of which for the frequency 
domains depend on the scattering angle 0, and hence on to and to'. Therefore, to be consistent, 
we apply Fourier decomposition to the redistribution matrix which contains the angle-dependent 



frequency domain information (see Anusha & Nagendra [111). This can be done as follows 



n 1 ^?? (m, x, x', 0, B) = ]T e ik * TZ^ Q KK (m, x, x' , 9, to', B), (11) 



k— — oo 

where the Fourier coefficients are given by 

n%^ K ' (m, x, x',9, to', B) = j* ^ c- ifc x TZ^ff' (m, x, x', 0, B). (12) 

The angle-dependent PRD functions i?n,m(x, x',0) are periodic functions of x with a period 2n 

because of which each element of the redistribution matrix TZqq,, is 27r-periodic. We remark that 
in the previous attempts on Fourier decomposition, the expansion of angle-dependent functions 
-Rn, 111(2;, x', 0) over (x~X') was traditionally used (see Domke & Hubeny [9(, Frisch [4(, Frisch [5J]). 
We show below that an expansion over \ of the angle-dependent redistribution matrix (as done 
in Anusha & Nagendra [ll|), provides a consistent way of including 'angle-dependent frequency 

domains' when performing angle-dependent PRD computations. The matrix elements TVqq,, are 
studied in detail in Sect. [3] Similar azimuthal Fourier expansions for the primary source term 
Gq (r) can be written as 

k— +00 

G%{t) = J2 ^ G Q K (^ ( 13 ) 

k— — 00 

with 

G Q (T) -\0 iffc^O. (14) 

Inserting the Fourier azimuthal expansions of TVqqh (m, x, x', 0, B) and Gq in Eq. (j9]), we obtain 
an expansion for SJ'q which can be expressed as 

fc=+oo 



k— — oo 



S^(r,x,to) = J2 ^ S}% K (r,x,e), (15) 



with 



sl k ]K (T,x,9)=G Q k)K (T)+ I dx> 



l,Q 

Q"=+K' 



+o ° , , r an' 1 



An ip(x) 
E E ^, K '{m,x,x',d,to',B)I Q ^\r,x',to'). (16) 



x 

A"=0,2 Q"=0 



Substituting from Eq. (IT51) for S\'q in formal solution we get 

k— +00 



i% K (T,x,n)= Y, eikx ! Q (T,X,0), 



(17) 



k= — 00 



where 



/^(r, x, 0) = / +3 ° e-( T '-^W/" S { Q k)K (r, x, 9) ^ dr', 



for /j, > 0, 



f'''-"s , am / -(t'-t) ¥ .(x)/m s(fe)K^ „ m ¥>(a:) j_, 



(r, x, 1 



^(r,x,, 



/'- 



dr , for fi < 0. 



(18) 



Thus from Eqs. ([T5j) and (|T7|) we get an expression for azimuthal Fourier source vector components 
as 



Mk)K 



S^(r,x,9)=G^(r) + 



_Mk)K, 



+00 



dx' 



dfi' 1 



47T y(x) 

y=+K' fe'=+oo 

;(k)KK' (rn „ „i a n' n\ ST J*'x' r&) K ' '1 



x E E Cf™. 1 ' 1 '^' 5 ) E e**%V*(T,x',0'). (19) 

K'=0,2 Q"=0 fe'=-oo 



Notice that the Fourier indices k and k! are not coupled to Q and Q" unlike in the case of decom- 
position over (x — x') ( see Frisch [4|, NS11). 

The advantage of working in real irreducible basis is that we can reduce the computational time 
by restricting the values of azimuthal Fourier index k to positive space using conjugate symmetry 
relations as shown below. This simplification is analytically complicated in complex basis. From 
Eq. (|T2|) we can see that the components TVqqii satisfy the symmetry relation 



n 



(k)KK' _ 
QQ" ~ 



K 



(-k)KK' 



Using the above relation in Eq. (fTTj) we get 



1£qqu (m, X, x', 0, B) — Re 



fc— +00 



J2 (2 - 6 k o)e ik ^TZ^ Q KK ' \m,x,x' ,$,& ,B) 



fc=0 



(20) 



(21) 



Notice that the Fourier series constitutes only the terms with fc ^ which is useful in practical 
computations. With this simplification and following as in Eqs. (IT51) - p^)) we get 



SJ'q (t, x, ft) = Re 



fc— +00 



]T (2-6 k0 )e^S { l k > K (T,x,9) 



L k=0 



and 



,K 



r i £(r,x,n)=Re 



'k—-\-oo 

J2 (2-4o)e ifex ^(r,x, ( 

. fc=0 



(22) 



(23) 



The §1 X (r, x, ft) now takes the form 



2"=+K' 



J oo J <P{ ) K , =02 Q"=0 

(24) 



xRe 



k'=+oo 

J2 (2~5 k , )e ik 'x'f Q k : )K '(T, X i,ei) 



k'=0 



3. Azimuthal Fourier components of the redistribution matrix elements 

In this section we present the azimuth angle dependence of the Fourier decomposed matrix 
elements of the redistribution matrix, namely TZ-qqi, (jn, x, x' , 9, Q' , B). From Eq. (|2"Tj) . it is clear 
that the value of k extends from to +oo. For numerical evaluation it is necessary to truncate this 
infinite series. Our studies show that the series can be truncated at k = 4. We compute Fourier 
components TZqq" numerically. This we do by numerical integration of TIqqh over the azimuth 
angle x using a Gauss-Legendre quadrature with 32 grid points between [0, 2ir]. 

Fig. [TJ shows the 7?-oo element of the redistribution matrix as a function of x' and \' for a given 
set of x, 9, 9' and the magnetic field parameters. The main feature is that the k = component 
is the dominant term. Even though k ^ terms depend sensitively on %', their magnitudes are 
several orders smaller than the k — component. For this reason, we can truncate the Fourier 
azimuthal expansion of the redistribution matrix elements to the fifth term itself without causing 
significant errors. In fact for practical computation one can truncate the series at k = 2 itself. This 
would help in rapid computation of the angle- dependent PRD problems in practical applications. 
However, in the theoretical studies presented in the paper we use k = 4. From Fig.[TJit follows that 
the higher order components show several harmonics as \' varies from to 27r. This behavior is 
confined to x 1 < 3 when x — 0. For larger values of x', the components approach the value zero. 
We have verified that the above conclusions remain valid for arbitrary choice of x, 9, 9' and the 
magnetic field parameters and for other combinations of K , K', Q, Q" . 

4. Scattering expansion method for Hanle effect with angle-dependent PRD 

In Sect. |3] we showed that the azimuthal Fourier expansion of the redistribution matrix (see 
Eq. (|21[) ) can be truncated to the fifth term. Thus, for k = 0, 1, 2, 3, 4, we obtain a finite set of 
54 coupled integral equations. The dimensionality of the problem increases to 54 complex quan- 
tities from 54 real quantities when we work in full space (i.e., —4 ^ k ^ +4). Thus working in 
positive half space is computationally advantageous. In this section we present an iterative method 
to solve this set of coupled equations. This method is based on Neumann series expansion of the 
components of the source vector contributing to the polarization. Sampoorna et al. [2j applied this 
method to solve the transfer problem with angle-dependent PRD in non-magnetic case and named 
it as SEM. These authors also show the efficiency of SEM over the core-wing-based polarized ap- 
proximate lambda iteration (ALI) method. NSI1 employed SEM to solve Hanle transfer problem 
with angle- dependent PRD. The results obtained by them showed slight inconsistency as compared 
to the results obtained from perturbation method (Nagendra et al. [3(). This might be due to the 
use of angle-averaged frequency domains (approximation III) to solve Hanle transfer problem with 
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Figure 1: Fourier azimuthal components of the type II redistribution matrix elements are shown as a function of x' 
and x' f° r x = 0, 8 = 7r/2, 9' = it and (r, 8b, Xb) = (1, 30°, 0°). The first five panels from left to right correspond 
to the real part of TZ 00 and the remaining panels refer to the imaginary part of 7?.qo 



angle-dependent PRD functions. Consistency in the results can be obtained by actually using the 
angle-dependent frequency domains (approximation II) . 

In SEM, first neglecting polarization we calculate Stokes /. We assume that Stokes / is cylin- 
drically symmetric and is given by the component Iq itself to an excellent approximation. This 
approximation yields k' = K' = Q" = in Eq. (fMf . The resulting component is the solution of a 
non-LTE unpolarized radiative transfer equation with the line source function given by 



^1,0 ( T ' x i' 



eB v 



+ OC 



cto' 



^~ 4t ^°o )0 °(m, x, x', 6, a', B) if °(r, x>, 9'). (25) 

47T <P(X) 



Eq. (|25p can be solved using a scalar ALI method based on a core-wing approach. Keeping only 
the contribution of J^ 0)0 on the RHS of Eq. ([21]) to the K = 2 coefficients and k = 0, 1, 2, 3, 4, 



(fc)2 



each component 5 1 ; q can be written as 



S$\t,x,9) 



(i) 



da;' 



dfi' 1 

47T <£>( 



J(fe)20 



(x) Q0 X ' ' 



0,tl',B)l(° )o (T,x',0'). (26) 



The superscript 1 stands for the single scattering approximation to the polarized component of 



the source vector. The corresponding radiation field 



'{k)2 

Q 



(i) 



for k — 0, 1, 2, 3, 4 is calculated 



with a formal solver and it serves as a starting solution for calculating the higher-order terms. The 
higher-order terms can be obtained by substituting for Ik,/ appearing in the RHS of Eq. (|24|). from 



f(k)2 



(i) 



We see that indices k, k' ', Q and Q' are now decoupled whereas they were coupled in the 



f(k)2 



(1) 



has 



case of Fourier decomposition over (x ~~ x')- Correspondingly the number of non-zero 

increased from 25 to 54 when changing from Fourier expansion over (% — x') to that over x- As a 
result the dimensionality of the problem has increased in the single-scattered solution computation. 
In the computation of higher order scattering terms, apart from keeping the coupling of (K = 
2, Q) components with other polarization components [K' — 2, Q"), we also keep the coupling of k 
with all other k' terms. We recall that in NS11 coupling of k with kl — terms were only retained. 
Thus S, X at order n are now given by 



'i.Q 



[S$tT t Z t 



(,,) 



J(fc)2 



(t,x,9) 



(1) 



+ Cx^ 



dx' 



dO' 1 



An (p(x) 



y=2 
Q"=0 



x E K%$(m,x,x?,9,n',B)Re{ ^ ( 2 -^'o)e 



k =+oo 



3 ife'x' 



f(k')2 



(r,x',e') 



(n-l) 



fc'=0 



.(27) 



From Fig. [T] it can be seen that k = component of the redistribution matrix elements dominate 
over the higher order terms (fc ^ 0). For this reason, despite a strong dependence of TZqq,, on 
x' , it is sufficient in the summation over kl to retain the leading term (namely k' = 0) in practical 
computations. The inclusion of higher order terms kl > do not affect the solutions significantly. 
Using only kl = term also saves great amount of computing effort. 

5. Results and discussions 

We compare the Stokes parameter / and the ratios Q/I and U / 1 computed from our present 
approach with the perturbation method of Nagendra et al. [3] and NS11 approach. The perturba- 
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Figure 2: Stokes profiles at /i = 0.112 computed using three numerical methods, namely the perturbation method 
(solid line), the NS11 approach (dotted line), and the self consistent approach used in this paper (dashed line). 
Model parameters are (T, a, e, r, B vo ) = (10, 10~ 3 , 10~ 3 , 0, 1). The magnetic field parameters are taken as (r, #3, 
Xb) = (lj 30°, 0°). The left panel shows the results computed using pure Ru function and the right panel is for a 
combination of Ru and Rm functions with Te/^r = 1. 
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tion method treats linear polarization as a perturbation to the scalar intensity, and computes the 
polarization in a two step process, wherein an accurately computed Stokes / is used as an input 
in evaluation of polarized scattering integral. In successive perturbations, the Stokes Q and U 
are computed more and more accurately until convergence is reached. NS11 approach is discussed 
in Sections Q] and 15.11 We consider self-emitting plane-parallel, isothermal atmospheres with no 
incident radiation at the boundaries. The slab models are characterized by (T, a, e, r, B VQ , Y^I^r 
), where T is the optical thickness of the slab, and Te is the elastic collision rate. The depolarizing 
collision rate D^ is set to T^/2. The plank function B„ is taken as unity at the line center. The 
polarizability factor W% is taken as unity (i.e, we consider aJ = 0— ► 1 — > scattering transition 
with J the total angular momentum quantum number). The vector magnetic field in the Hanlc 
scattering problem is defined through the field strength parameter T = gco^f^R, with g the Landc 
factor of the upper level and «£ the Larmor frequency; the field inclination (Ob, Xb) defined with 
respect to the atmospheric normal. For angle and frequency discretization we have used quadratures 
of the same order as those used by NS11. Therefore we do not elaborate here on the computational 
aspects. 

5.1. A comparison with previous approaches to solve the angle- dependent Hanle transfer problem 

In this section we present the Stokes profiles of the lines formed in a magnetized slab scattering 
according to Hanlc PRD matrix formulated by Bommier [lCf . In the so called approximation II and 
III of Bommier [10( , the switch over from the Hanle phase matrices (in the core) to the Rayleigh 
phase matrix (in the wings) is achieved through the use of angle-dependent and angle-averaged 
frequency domains respectively. It is shown by Bommier }10| that the use of frequency domains 
simplifies the numerical evaluation of the redistribution matrices. 

It is natural that in AD radiative transfer computations involving AD functions, one should use 
Approximation II involving AD functions. However, as already discussed in Sects. Q] and [4] in NS11 
AA domains were used in computing the AD redistribution matrix. Although such an approach is 
inconsistent, it provides a rapid means of solving the Hanle AD PRD problems. In the present paper 
we test their approach by actually using AD frequency domains while computing AD redistribution 
matrix (which is fully consistent). In Fig. [21 we show a comparison of results obtained by the NS11 
and the present approach, along with those obtained from the simple perturbation method of Na- 
gendra et al. [3J], which is also consistent, like the present approach. The left panels in Fig. [2] show 
the Stokes profiles computed using pure R\\ function. One can clearly see that the NS11 approach 
differs from the present approach particularly in the frequency range 3 < x < 5. The present 
approach and perturbation method give same results. The impact of the approximation used in 
the NS11 is more severe on the Stokes U parameter. The right panels in Fig. [5] show the results 
computed using the same model as in the left panels, but for the introduction of elastic collisions (a 
combination of Ru and Rm). One can clearly see that the differences between NS11 and the present 
approach still exist, although the collisions decrease these differences. The present approach, unlike 
the approximate treatment followed in NS11, thus provides a self-consistent approach to compute 
the redistribution matrix (i.e., the use of AD domains to compute AD redistribution matrices), 
at the same time requiring manageable computing resources . This has practical implications in 
realistic modeling of the observed Stokes profiles. 

Table [1] shows a comparison of computing resources required by three numerical methods. Com- 
pared to the perturbation method that requires large computing time (112 minutes to obtain a 
solution), the approximate method of NS11, and the present method are less expensive. In spite of 
being inconsistent, the approximate method of NS11 requires far less computing memory compared 
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Figure 3: A comparison of emergent Stokes profiles computed by three numerical methods discussed in the text. The 
profiles are presented for fi = 0.112 and for a vertical magnetic field (r, 9g, \b) = (1> 0°, 0°). Different line types 
represent same cases as in Fig. [2] The slab model parameters are (T, a, e, r, B VQ , Te/^r) = (2 x 10 4 , 10 -3 , 10 — 3 , 
0, 1, 1). 
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Table 1: Comparison of CPU time taken by different methods for radiative transfer computations. The model 
parameters used for the computation are (T, a, e, r, B VQ , Ve/Tr) = (2 x 10 6 , 10~ 3 , 10~ 3 , 0, 1, 1) 





Time (minutes) 


Memory 


Present approach 

Perturbation method 

NS11 approach 


28 
112 
33 


13GB 

7.6GB 
226MB 



to the other two methods. This is because in NS11 the polarized transfer equation is solved in 
a azimuth independent Fourier basis, which thereby avoids introducing azimuth angle grids. The 
present approach requires larger memory because we now need to discretize the azimuth angle \' an d 
store the huge matrix 'R-nn" ( m i %> %' ■> $> $~2', B). The memory requirement of the present approach 
is even larger than that of perturbation method because the former involves solving 54 coupled in- 
tegral equations in Fourier basis, while the later involves solving 3 coupled integral equations in 
Stokes basis. However, unlike the present approach the convergence is not always guaranteed in the 
perturbation method, as angle-frequency coupling is more intricate in Stokes basis than in Fourier 
basis. 

5.2. The vertical field Hanle effect 

It is expected that when the magnetic field is parallel to the symmetry axis of the slab (the 
atmospheric normal) , the Hanle effect should vanish. In other words the Stokes U parameter should 
be zero in this case. This characteristic behavior is satisfied when we work with angle-averaged 
redistribution functions. When angle- dependent redistribution function is used this behavior is not 
satisfied. In other words, Stokes U does not vanish, in spite of the field being vertical, as long as the 
AD redistribution function is used. The non-zero emergent Stokes U is due to coupling of Stokes U 
to Stokes / through the components TIqq, for k =/= 0. The reason for this unexpected behavior is 
also discussed by Frisch et al. [18| and numerically demonstrated in Nagendra et al. |3j . We revisit 
this interesting problem in Figs. [3] and 01 In Fig. we show the Stokes profiles computed using 
NS11 and the present approach, and compare it with the results from perturbation method. The 
differences between the NS11 and the present approach are drastic in Stokes U parameter. The 
Stokes U in the frequency range 3 < x < 5 has opposite signs. The present approach produces 
Stokes U consistent with the perturbation method. The NS11 approach for approximation II seems 
to be inadequate in computing U in this particular problem. Such a large difference between the 
results obtained from NS11 approach and the perturbation method prompts us to conclude that it 
is safer to use AD frequency domains to solve AD Hanle transfer problems. 

Fig. 2] shows the center to limb variation of linear polarization for a vertical magnetic field. 
The intensity exhibits the characteristic limb darkening in the line core and limb brightening in 
the wings. The Q/I shows limb brightening throughout the line profile. The dependence on /i is 
non-monotonic in the core region in U/I. Since the angle- dependent functions become azimuthally 
symmetric at the disk center the U/I approaches zero as the line of sight approaches the disk center. 
In the line wings, U/I tend to zero for all /i's which is due to the Rayleigh scattering in the line 
wings, that produces U/I = by axisymmetry. 

5.3. The Hanle effect with micro-turbulent magnetic field 

It is known that the presence of a weak turbulent magnetic field in the solar atmosphere can 
be detected using Hanle effect. In the case of CRD, Frisch et al. [8] showed that the polarization 
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Figure 4: Stokes 7 and ratios Q/I and {//J for different emergent angles (9 = cos -1 fi) computed for a vertical 
magnetic field. The atmospheric model is same as in Fig. [3] Different line types are: solid line - fi = 0.025, dotted 
line - fj, = 0.129, dashed line - /i = 0.297, dot-dashed line - fJ. = 0.50, dash-triple-dotted line - fi =0.702, long-dashed 
line - fi = 0.871, and thick solid line - fi = 0.974. 
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Figure 5: Stokes Q/I for emergent angle fi = 0.112 computed for the case of micro-turbulent magnetic field. The 



atmospheric model is (T, a, e, r, B VQ , T~e/^r) = (10, 10 



io- 



0, 1, 0). Panel (a) shows the comparison of the 



present approach (solid line) and NS11 approach (dotted line). Panel (b) shows the comparison of angle-averaged 
and angle-dependent cases. Different line types are : non-magnetic case - solid lines; deterministic magnetic field - 
dashed lines; and micro-turbulent magnetic field - dot-dashed lines. The thick and thin lines represent angle-averaged 
and angle-dependent cases respectively. For the micro-turbulent magnetic field case, T = 1 and for the deterministic 
magnetic field case, (r, 6b, Xb) = (1, 30°, 0°). 
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obtained using Hanle effect is quite sensitive to the choice of field strength distribution and in 
general, micro-turbulence is a safe approximation to represent weak turbulent magnetic fields. For 
a micro-turbulent magnetic field, the scale of variation of the field is small compared to the mean 
free path of the photons, and this allows to replace all the field dependent physical parameters by 
their averages over the magnetic field vector probability density function (PDF). In our problem this 
condition leads to the averaging of the magnetic kernel Mqq, (m, B) over the magnetic field vector 
PDF. In this paper we use a PDF corresponding to the isotropic distribution of field orientation {6b, 
Xb) and a single value of the field strength. As shown in Stenflo [19j, [20j the Hanle problem with 
this choice of PDF reduces then to a resonance polarization problem, with a modified value of Q/I . 
In other words, the micro-turbulent averaged magnetic kernel namely < Mqq, (m, B) > becomes 
diagonal, and only < A/Jo (m, B) > element is of relevance. The explicit form of < A/Jq (m, B) > is 
given by (see Appendix B of Frisch et al. [8|) 



rtmBW ., 2 T 2 (m) 2 AT\m) 



<M^{m,B) >=!-- - ■ -y - -- - ■ -A-J (28) 



where Y(m) denotes the Hanle T-parameter in different frequency domains. We refer the reader 
to Eqs. (89) of Bommier [10( for the explicit form of T(m) in different frequency domains where 
Hanle effect is operative. In frequency domains where Rayleigh scattering is present, T(m) = 0. 
Fig. 0(a) shows the ratio Q/I obtained using the present approach and the NS11 approach in the 
presence of micro-turbulent magnetic field. The difference between the results obtained using these 
different methods mainly exists in the transition region 3 < x < 5. Fig.[5](b) shows the comparison 
of Q/I profiles computed using angle-averaged and angle- dependent redistribution matrices. For 
AA computations, we use redistribution matrices computed in AA domains (approximation III 
of Bommier [l0|), and for AD computations the corresponding AD redistribution matrices are 
computed in AD domains (approximation II of Bommier [10]). Three different sets of results are 
presented in Fig. [5](b), namely non- magnetic, deterministic field, and micro-turbulent field results. 
We see that in the case of Hanle effect with micro-turbulent magnetic field, there is depolarization in 
the line core and in the near wing frequencies as compared to the other two cases. Our studies show 
that in the case of Hanle effect with micro-turbulent magnetic field, the differences between angle- 
averaged and angle-dependent results are prominent in thin slab cases and reduce considerably 
in thick slab cases. These differences between AA and AD results, especially in the 3 < x < 5 
region (apart from the line core) , were already noticed by Nagendra et al. [3j , for the deterministic 
magnetic field case. It is interesting to note that such differences get enhanced in the presence of a 
micro-turbulent magnetic field. This is probably because of the localization of line photons within 
the micro-turbulent scattering eddies, resulting in a relatively larger number of scattering. These 
differences are present both in the CRD scattering mechanism (see Frisch et al. [8(), as well as PRD 
as can be seen from Fig. [5] (b) . 

6. Conclusions 

We have solved the angle- dependent Hanle scattering problem using angle-dependent partial 
frequency redistribution (PRD) theory (approximation II of Bommier [10(). This computationally 
expensive problem is solved using an iterative method based on the Neumann series expansion 
(SEM). Following Anusha & Nagendra [16j we decompose the Stokes parameters in terms of az- 
imuthally symmetric Fourier coefficients, by expanding the Hanle redistribution matrix in terms of 
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the radiation azimuth \. Only such a decomposition allows the use of angle-dependent frequency 
domains for solving angle-dependent Hanle scattering problems. In contrast a decomposition, based 
on the expansion in terms of (x ~ x')> as done in NS11, does not allow the use of angle-dependent 
domains to solve angle- dependent Hanle scattering problem. For this reason a simpler approach 
was suggested by NS11 which used angle-averaged frequency domains (approximation III of Bom- 
mier [10]), to solve angle- dependent Hanle transfer problem. We show that their approach does 
not always hold good. We have carried out a numerical study to show the differences between 
the solutions obtained by NS11 approach and the self-consistent approach used now in this paper. 
The U/I profiles in particular show significant differences in the core to wing transition region 
(3 ;$ x < 5) of the line. The special case of vertical field Hanle effect is considered as a case study, 
and the differences between the NS11 and the present approach are examined. It is shown that 
the present method offers a self-consistent and accurate method of solving the difficult problem of 
AD partial redistribution with Hanle scattering. The interesting behavior of Q/I profiles in the 
presence of micro-turbulent magnetic fields is also examined. We show that the differences between 
angle-averaged and angle-dependent solutions are enhanced by the presence of a micro-turbulent 
field. The differences are noticed in both the line core and near wing regions (3 < x < 5). 
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